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I. INTRODUCTION 



In the past few years there has been considerable work on non- Hermit ian Hamiltonians. 
Among this kind of Hamiltonians, much attention has been focused on the investigation of 
properties of so-called 'PT-symmetric Hamiltonians. Following the early studies of Bender 
et al. [1], the VT— symmetry formulation has been successfuly utilized by many authors 
[2-8]. The PT-symmetric but non-Hermitian potentials could have real spectra even it is 
non-Hermitian. Non-Hermitian Hamiltonians with real or complex spectra have also been 
analyzed by using different methods [3-6,9]. Non-Hermitian but PT-symmetric models 
have applications in different fields, such as nuclear physics [10], condensed matter [11] and 
population biology [12]. 

When a particle is in a strong potential field, the relativistic effect must be considered, 
which gives the correction for non- relativistic quantum mechanics [13]. Taking the rela- 
tivistic effects into account, a particle in a potential field should be described with the 
Klein-Gordon (KG) and Dirac equations. In recent years, there has been an increased inter- 
est in finding exact solutions to Schrodinger, K-G, Dirac and Salpeter equations for various 
potential schemes [13-27]. The problems that can be exactly solved for the KG and/or 
Dirac equations are seldom except a few examples, such as hydrogen atom and electrons 
in a uniform magnetic field. Recently some authors solved such relativistic equations for 
some potentials. §im§ek and Egrifes [15] have presented the bound-state solutions of the 
one-dimensional (ID) Klein-Gordon (KG) equation for PT-symmetric potentials with real 
and complex generalized Hulthen potential. Moreover, Egrifes and Sever [16] investigated 
the bound state solutions of the ID Dirac equation with PT-symmetric real and complex 
forms of generalized Hulthen potential. Yi et al. [19] obtained the energy equations in the 
KG theory with equally mixed vector and scalar Rosen-Morse-type potentials. In recent 
works, we have solved the spinless ID Salpeter equation analytically for its exact bound 
state spectra and wavefunctions with real and complex forms of the PT-symmetric general- 
ized Hulthen potential [20]. We have also investigated the bound state solutions of the ID 
KG equation with real and complex forms of the generalized Woods-Saxon (WS) potential 
[21]. 

In the present work, we present a new procedure to construct solution to the (1+1) dimen- 
sional Dirac equation with gauge invariant (minimal) vector coupling. Using the Nikiforov- 
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Uvarov (NU) method [28], we are set to obtain the bound states solutions (relativistic energy 
spectrum and the two-component spinor wave functions) for particles trapped in a spheri- 
cally symmetric, generalized Woods-Saxon potential [21,22] which possesses a "PT-symmetry 
as well. Changing the values of the potential parameters from real to pure imaginary, we 
obtain Hamiltonians that may or may not be PT-symmetric. Further, by making the same 
parameter change in the energy spectrum and two-components spinor wave functions, we 
obtain solutions for the new trigonometric and periodic potential (real or complex) forms. 
The corresponding energy spectra are either real or complex. The Dirac problem with a 
"PT-symmetric or a non-PT-symmetric imaginary generalized Woods-Saxon (WS) poten- 
tial is mapped into the exactly solvable problem in which one may apply the NU method to 
generate possible real energy spectra [28] . 

The paper is structured as follows: In Section [Til we briefly introduce the basic concepts 
of the hypergeometric method (NU). Section [TTTJ is devoted to the solution of the Dirac 
problem to obtain the exact bound state energy spectrum for real and complex cases of 
generalized WS potentials and the lower and upper spinor components eigenfunction by 
applying the NU method. The results of this relativistic study are discussed in Section ??. 

II. THE HYPERGEOMETRIC METHOD 

The basic equations of the theoretical background are similar to those given in [15,16,20- 
25]. According to a breif description of this method [28], the Dirac equation can be trans- 
formed to the following generalized equation of hypergeometric type after employing an 
appropriate coordinate transformation, s = s(r), 

where a(s) and a(s) are polynomials, at most of second- degree, and f(s) is a first-degree 
polynomial. In order to find a particular solution to Eq. (1), we use the following wave 
function 

4>n(s) = <p n (s)y n (s). (2) 

This reduces Eq. (1) to an equation of a hypergeometric type 
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a(s)y:(s) + r(s)y' n (s) + \y n (s)=0, 
which demands that the following conditions be satisfied: 



(3) 



°m=«'Wy (4) 

r(s) = r(s) + 2tt(s), t'(s) < 0, (5) 

and A is thus a new eigenvalue equation for the second-order differential equation becomes 

A = A n = -nr'(s) - " ^ ~ l) a"(s), n = 0, 1, 2, .... (6) 

The polynomial r(s) with the parameter s and prime factors show the differentials at first 
degree be negative. It is worthwhile to note that A or A n are obtained from a particular 
solution of the form y(s) = y n (s) which is a polynomial of degree n. The other part y n (s) 
of the wavefunction (2) is the hypergeometric-type function whose polynomial solutions are 
given by Rodrigues relation 

Vn{s) = W)^ [an{S)p{S)h (7) 
where B n is the normalization constant with the weight function p(s) must satisfy the 

following condition 

w\s) - (^j w(s) = 0, w(s) = a(s)p(s). (8) 

On the other hand, in order to find the eigenfunctions, (f> n (s) and y n {s) in Eqs. (4) and (7) 
and eigenvalues A n in Eq. (6), we need to calculate the functions: 



cr'(s) — t(s)\ 2 _ . , , . , 



A = A; + 7r'(s). (10) 

In principle, since 7r(s) has to be a polynomial of degree at most one, the expression under 
the square root sign in Eq. (9) can be arranged to be the square of a polynomial of first 
degree [28]. This is possible only if its discriminant is zero. Thus, the equation for k obtained 
from the solution of Eq. (9) can be further substituted in Eq. (10). In addition, the energy 
eigenvalues are obtained from Eqs. (6) and (10). 



III. SOLUTIONS OF THE QUANTUM SYSTEM 



The exact solution of the relativistic wave equations is of much concern in quantum 
mechanics. Many works appeared in the recent years in that direction towards obtaining 
exact solution of some relativistic wave equations for certain potentials of physical interest 
(cf. [15,16,20-22] and references therein). For a spinless particle of rest mass m and total 
energy E n i, the ID time- independent Dirac equation with any given interaction potential 
V(x) in the vector coupling scheme (choosing the natural atomic units H — c — 1) is 
[14,16,29-31] 




+ [E n - V{x)\ 



1 

1 




V>„ g (x) = 0, (11) 



Zi n (X) 

where the spinor radial wave function ^ nq (x) = j ] is normalized with u nq (x) 

w nq (x) 

(w nq (x)), the upper (lower) spinor component (the x range is from — oo to oo on the full-line 
problem). The interaction among nuclei is commonly described by using a potential which 
consists of the Coulomb and the nuclear potentials. It is usually taken in the form of Woods- 
Saxon (WS) potential. Here we take the following Hermitian real-valued ID generalized WS 
potential which is specified by the shape (deformation) parameter, q, [22,32,33] 

V q (x) = -V —^ -, a = 1/a, r-Ro^x, q > 0, #0 > a, (12) 

1 + qe ax 

where r refers to the center-of-mass distance between the projectile and the target nuclei 
(the r range is from to oo). The relevant parameters of the inter-nuclear potential are 
given as follows: R = r^A 1 ^ is to define the confinement barrier position value of the 
corresponding spherical nucleus or the width of the potential, A is the target mass number, 
7o is the radius parameter, the field strength Vq controls the barrier height of the Coulombic 
part, a is the surface diffuseness parameter has to control its slope, which is usually adjusted 
to the experimental values of ionization energies. Note further, q is a shape (deformation) 
parameter, the strength of the exponential part other than unity, set to determine the shape 
of potential and is arbitrarily taken to be a real constant within the potential. It should be 
noted that the spatial coordinates in the potential are not deformed and thus the potential 
still remains spherical. 
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It is worth to state that under radial coordinate transformation, r — > r + A, then the 
generalized WS potential in Eq. (12) changes into the standard WS potential (q = 1) but 
with the displacement parameter A satisfies the expression exp(A/a) = q and with a field 
strength Vq = Vq exp(— A/a) [32]. The sense of generalization or deformation of the potential 
becomes clear. For completeness, it could be stated that if A is positive (corresponding to 
q > 1) then one may need to impose the condition on the choice of A, that is, |A| <C R Q . 

Obviously, for some specific q values this potential reduces to the well-known types, such 
as for q = to the exponential potential and for q = — 1 and a = d^ 1 to the generalized 
Hulthen potential (cf. [20,34] and the references therein). 

For the given generalized WS potential, Eq. (11) decomposes into the following two 
components: 

/ d e ~ ax \ 

\dx + Enq + V ° 1 + qe~ ax J Unq ^ = mw "<?( x )> ( 13 ) 
/ d e ~ ax \ 

\dx + Enq + V ° l + qe~ ax ) Wnq ^ = mu m( x )- ( 14 ) 
It should also be noted that the coupled differential equations allow finite bound states 
(real) solutions as u nq (r- — > oo) — > and w nq (r — > oo) — > 0. The above set of coupled 
equations can be reduced to a second-order differential equation. Hence, combining the last 
two equations, this provides the following second-order Schrodinger-type equation for upper 
(lower) components, respectively, as 



O) + 



Kq( X ) + 



e -2ax 



(1 + qe ax ) 2 1 + qe 

p —2ax p—oi: 
En q + V* - ■ —— 2 +V 2 



1 + qe~ ax ) 2 1 + qe 



Ung(x) = 0, (15) 
w nq (x) = 0, (16) 



with the following definitions 



E nq = E 2 nq - m 2 , V 1 = V 2 - tqV /a, V 2 = 2E nq V + iV /a, (17) 

and V* (y 2 *) is the complex conjugate of V\ (V2), respectively. 

Using the NU method, we are set to obtain bound states solutions (relativistic en- 
ergy spectrum and spinor wavefunctions) of a spin-zero particle for a three parameter 
{V 0: q, ajgeneralized WS potential. We employ the following dimensionless transformation 
parameter, s(x) = (1 + qe' 01 *)' 1 , which maintains the transformed wavefunctions finite, on 
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the boundary conditions (i.e., < r < oo — > — oo <:r<oo^0<s<l) [21,22]. Hence, 
Eq. (15) is reduced into the generalized equation of hypergeometric type which is given by 
Eq. (1): 



d 2 u nn {s) . 1 - 2s du nn (s) . a 2 . , , 

u nq {s) = 0, (18) 

Q Q 



^ + 5(l-s) 2 + ^(l-s) 



(is 2 (s — s 2 ) rfs (s — s 2 ) 2 
where we have set u nq {x) = u nq (s). Therefore, with the dimensionless definitions 

-e 2 = a 2 E nq > 0, (3 2 = a 2 V 2 /q, 7 2 = a^/g 2 , (E 2 nq < m 2 , (3 2 > 0, 7 2 > 0), (19) 

for bound states (i.e., real e 2 ), one can arrive at the simple hypergeometric equation given 
by 

Before further proceeding, it is necessary to compare the last equation with Eq. (1) to obtain 
the following polynomials: 

r(s) = 1 - 2s, (t(s) = (s - s 2 ), (t(s) = s 2 7 2 - s ((3 2 + 2 7 2 ) + (3 2 + 7 2 - e 2 . (21) 
The substitution of the above expressions into Eq. (9), together with cr'(s) = 1 — 2s, gives 



tt(s) = ±i y/s 2 ( 7 2 + fc) - s 03 2 + 2 7 2 + k) + (3 2 + 7 2 - e 2 . (22) 

It is taken into consideration that the discriminant of the square root sign has to be zero. 
Hence, the second-order equation under the expected roots are obtained as fc lj2 = (3 2 — 2e 2 ± 
2eb, where b = y/e 2 — (3 2 — 7 2 = — (n + 1 + iaV /q) — e with n = 0, 1, 2, • • • . In this case, 
substituting these values for each k into Eq. (22), the possible solutions are obtained for 
7r(s) are: 

{(b-e)s-b; for kl = p 2 - 2e 2 + 2eb, 
tt(s) = ± < (23) 

I {b + e) s - b; for k 2 = f3 2 - 2e 2 - 2eb. 
For bound state solutions, it is necessary to choose one of the four possible forms in the last 
equation. Therefore, the most suitable form is established by 
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7r(s) = -(b + e)s + b, k = (3 2 - 2e 2 - 2e6. (24) 

The following track in this selection is to achieve the condition t\s) < in Eq. (5), which 
can be obtained as 

r(s) = -2(1 + b + e)s + 1 + 26, 

t'(s) = -2(1 + b + e) = 2(n + iaVb/g). (25) 

A particular solution can be calculated by using Eqs. (6) and (10). Consequently, this 
solution is obtained as 

A = - 7 2 - (b + e)(b + e+l), 



A n = n 2 + n + 2n (e + b) . (26) 
After setting A n = A and solving for E nq , we find the Dirac exact energy spectra as 

V 



E. 



nq 



2q 



± Kn(q, a, V 



\V 2 + Kl(q,a,V ) ^ 



K n (q, a, V ) = - [iV + qa{n + 1)] , n = 0, 1, 2, • • • (27) 

We should point out that Egrifes and Sever [16] have recently obtained a similar expres- 
sion to Eq. (27) for the case of the Hulthen potential. Let us now find the corresponding 
wavef unctions. We have seen that the energy expression has a complex form for the poten- 
tial under study. We look for the complex generalized WS potential forms that have real 
spectrum. 

Moreover, the restriction which gives the critical coupling value leads to the result 

1 



n < 



qa 



4q 2 m 2 _ yi _ iVo _ 1 



(28) 



that is, there are only finitely many eigenvalues. In order that at least one level might exist, 
it is necessary that the inequality 



qa + iV ^ \jAq 2 m 2 - V 2 , 



(29) 
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is fulfilled. As can be seen from Eq. (28), there are only two lower- lying states for the Dirac 
particle of mass unity when the parameters a — 1, q — ±1 for any given Vq : 



n 



<±(y/4- V 2 - iV^j - 1. (30) 



(i) Choosing q — 1, the potential form (12) is reduced to the shifted WS potential: 

nx) = -Vo + T ^, (31) 
and then it's energy spectra yield 

E n = -— ± [iV + a(n + 1)1 , 

2 ;j ^ 2 + [^ + a(n + l)] 2 4 

(ii) Choosing g = —1, the potential form (12) is reduced to the shifted Hulthen potential: 

= Vo~ (33) 
and then the resulting energy eigenvalues become 



171 ^ ^ 

2 - 7 , n = 0,l,2,... (32) 



E n = ^±\iV -a(n + l)] 



171? ! 

„ = 0,1,2,... (34) 



2 /J ^\/ 2 + [^ -a(n + l)] 2 4 

(iii) For the case g — ► 0, the potential expression (12) is reduced to the exponential potential: 

V(x) = -V e~ ax , (35) 

the eigenvalues expression (27) does not give an explicit form, i.e., the NU method is not 
applicable to the exponential potential (34). 

Note that for this potential there is no explicit form of the energy expression of bound 
states for Schrodinger [18], KG [35] and also Dirac [17] equations. 

In addition, it can be seen easily that while the field strength Vo — > 0, the energy eigen- 
values yield: 

E n = ±l^m 2 - (n + l) 2 a 2 , n = 0, 1, 2, • • • (36) 

Note that in the above equation there exist bound states for the ground and excited states 
(n = 0, 1) which are E = ±y / 3m/2 and E\ = 0, respectively, for positive q values and 
a = A c , where \ c = 1/m denotes the Compton wavelength of the Dirac particle. Otherwise, 
there are no bound states for n > 2 states. 
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On the other hand, for the same value of a and negative q values when Vq — > 0, all energy 
eigenvalues go to zero. If the value of q is increasing, all positive bound states go to zero, 
from (27), asymptotically. 

An inspection of the energy expression given by Eq. (27), for any given a, shows that we 
deal with a family of generalized WS potentials. The sign of V does not effect the bound 
states. The spectrum consists of complex eigenvalues depending on q. As we shall see the 
role played by the range parameter a is very crucial in this regard. Of course, it is clear 
that by imposing appropriate changes in the parameters {a, V , q} , the energy spectrum in 
Eq. (27) for any modified parameter can be also calculated by resolving Dirac equation for 
every parameter change. 

Let us calculate the wavef unctions. Inserting, n(s) and a(s) in Eq. (4) and consequently 
solving the resulting first-order differential equation, we find 

Ms) = s b (l - s) e . (37) 

In addition, to find the function, y nq (s), which is the polynomial solution of hypergeometric- 
type equation, we multiply Eq. (3) by p(s) so that it can be written in self-adjoint form 
[28] 

(<j(s)p(s)y' nq (s)y + \p(s)y nq (s) = 0, (38) 
where p(s) satisfies Eq. (8), which gives 

p(s) = s 2b (l-sr. (39) 
The second eigenfunction can be obtained by Eq. (7) as 



d 



n 



y m {s) = D nq s- 2b {l - S )" 2 ^ [s n+2b (1 - 8 y*"\ , (40) 

where D nq is a normalization constant. In the limit q — > 1, the polynomial solutions of 
y n {s) are expressed in terms of Jacobi Polynomials, which is one of the classical orthogonal 
polynomials, with weight function given by Eq. (39) for s E [0, 1] , giving y n (s) ~ P^ b ' 2t \\ — 
2s). Obviously, the radial wave function u ng (s) for the s-wave can be obtained by substituting 
Eqs. (37) and (40) into Eq. (2) as 
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u n (s) = N nq s\l - S yPi 2b *\l - 2s), (41) 

where s(r) = (1 + e -(»--Rb)/«)-i > R > a and N nq is a new normalization constant. Further, 
by using the differential and recursion properties of the Jacobi polynomials [36], the lower 
spinor component can be obtained from Eq. (13) as 



mw ng (s) = N nq s b (l - s) e \ E nq - iae + 



la In + 1 H (1 — s) 

q \ aq 



| P (^) (1 _ 2s) 



+ N ng ia (n + l + ^ ) s\l - S yP^ b+1 ^ +1 \l - 2s). 
V qa J 



(42) 



Notice the well behavior of the wave function at infinity. As an example, the ground state 
wave function behaves like 

u (s -> 0) -> 0, m (s -> 1) -> 0, (43) 

and thus representing a truly bound-state solution. Further, we make use of the fact that 
the Jacobi polynomials can be explicitly written in two different ways 



pM {z) = 2 -n ±(-1)^ ^ + P\( n + A (1 _ Z) n-V (1 + z)P (44) 

p=o \ P / \ n V) 



T(n + p + l) " /n\ r(n + p+z/ + r + l) ^ - 1 
n!r(n + p + z/+l)^ Vr/ T(r + p+l) 



(45) 



where (™) 



r(r+i)r(t t 1 -7-+i) • Using Eqs. (44) and (45), we obtain the following two 



r!(n— r)\ 

explicit expressions: 



Pi 2b ' 2e \l - 2s) = (-l) n r(n + 2b + l)r(n + 2e + 1) 



x 



" (-1)V P 

- P) !r G° + 2e + l)r(n + 26 - p + 1) 



s™- p (l - s) p , 



P^(l-2s) = 



T(n + 26+ 1) 
r(n + 2b + 2e + 1) 



1)5 



r=0 



-l) r gT(ra + 26 + 2e + r + l) 
r!(n-r)!r(26 + r + 1) 



s' . 



(46) 
(47) 



1 = N 2 nq {-1) 



; r(n + 2e + l)r(n + 26 + l) 



2 n 



E 



(_l)P g n-p 



r(n + 2e + 2b + 1) 1 ^ p\(n - p)!r(p + 2e + l)T(n + 2b - p + 1) 
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X 



(-l)Vr(n + 2e + 26 + r + l) l 
J- r!(n-r)!r(2fo + r + l) /^^'^' ^ 



where 



r) = / s «+ 2 ^-^(l - (49) 
o 

Using the following integral representation of the hypergeometric function [36] 

}s"°-\l - 8)-»-°°-\l- q 8)-*>d8 = 2 Fl (a ,(3 : ™ q ) T{a0) l^- a0) , 

[Re( 7o ) > Re(a ) > 0, |arg(l - q)\ < tt], (50) 

which gives 

i 

2 F 1 (a ,[3 :a + l;q)/a = J s^'^l - qs^ds, (51) 

o 

TP/ a \ r(7o)r(7o - a -Pq) 

2 F 1 (a ,(3 : 7o; q) = ^ r— ; 

T(7o - a )r(7o - Po) 

[Re( 7o - a - (3 ) > 0, Re( 7o ) > Re(/3 ) > 0, (52) 
for q = 1. Setting ao — n + 2b + r — p + 1, /3q — —p — 2e, and 70 = «o + 1, one gets 

T ( x _ 2^1 (ao, Po ■ 7o; q) _ (n + 2b + r - p + l)\{p + 2e)\ 

In ^ T) a (n + 2b + r-p+l)(n + 2e + r + 2b + l)V { ' 

In view of the above complex energy spectra (27), it will be of interest to see how com- 
plex potential form of Eq. (12) would effect this result. Therefore, we shall change the 
values of the potential parameters (a, q, Vq) from real to pure imaginary (complex) to ob- 
tain Hamiltonians that may or may not be VT— symmetric. Hence, we also make the same 
parameter change in the energy spectra (27) and the upper and lower spinor components 
of the wavefunction Eqs. (41) and (42) respectively. The resulting non-Hermitian complex 
potential could have real energy spectra. To this end we consider the complexified forms of 
the generalized WS potential. 
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A. Non-Hermitian VT — symmetric new trigonometric and periodic potential form 



Let us consider the case where at least one of the potential parameters be complex. In 
this case, a is taken to be a complex parameter (i.e., a — > ia). Consequently, the potential 
in Eq. (12) transforms into the form 

Vq ^ = ~ q 2 + 2qll(ax) + l ^ + °° S ^ ~ lsin ^ = K(~ x )> ( 54 ) 
which is a VT — symmetric but non-Hermitian. We note that the transformed potential in Eq. 
(54) has a trigonometric and periodic form. Obviously, the last form of potential forms have 
no physical relation what so ever with the WS potential in (12), any of it's generalizations, 
it's well-known behavior or properties. As we have seen, simple mathemetical manipulations 
have not to cloud our intutive judgement and conceptual physical understanding. From 
a mathematical prospective, it might be possible that one can use calculus to think of 
the hyperbolic or exponential functions as another form of trigonometric ones. However, 
physically these functions, if considered as potentials, are dramatically different. We should 
not refer to either potential in Eq. (54) as WS-type. Nevertheless, this new complex potential 
embodies their periodic, trigonometric and VT — symmetric behaviors. Hence, this type of 
potentials (54) has real spectrum given by 



E nq = -^-±(V + q(n + l)/a) 



m? 



\2 ' 



(55) 



2? V " '\*q 2 V 2 -(V + q(n+l)/a) 

if the following restriction Aq 2 m 2 < V 2 — (Vo + q(n + 1) /a) 2 being achieved. The critical 
coupling value is 

q(n + 1) 2qam 2 

which leads to the following condition: 



V + qa- JV§ - 4q 2 m 2 V + qa + a/V 2 - 4q 2 m 2 

< n < , 

qa qa 

meaning that the number of real eigenvalues are finite. Further, the corresponding radial 
wave function u nq {s) for the s-wave could be determined as 

u ng (s) = N nq s c (l - S )-Pf^)(l - 2s), (56) 
13 



mw nq (s) = N nq s c (l - s) ie I E nq + iae + 



h a [n + H (1 - s) 

q V aq 



|pi 2c ' 2ie )(l-2s) 



N nq a ( n + l + ^\ s c (l - s )^p(Wie + i) {1 _ 2s)) 



(57) 



where c = —(n + 1 + aK)/g) — ie, and s(r) = (1 + qe~ l ^ r ~ Ro ^ a )~ 1 . 

For the sake of comparing the relativistic and non-relativistic binding energies, we need 
to solve the ID Schrodinger equation for the complex form of the generalized WS potential 
given by Eq. (54). Employ a convenient transformation given by s(r) = (l + ge-^-* )/ )- 1 , 
0<r<oo^0< s(r) < 1, we obtain [21] 



1 - 2s 



-f3 2 s + /3 2 - 6 21 
(s - s 2 f 



^ng(-s) = 0, 



(58) 



for which 



f(s) = l-2s, a(s) = s-s 2 , a(s) = -(3 2 s + (3 2 - e 2 , 



9 2ma ,_, 
e 2 = 



_ 2 ^ng 



2ma z 
~ti*q~ 



V (E nq <0, (3 2 > 0). 



(59) 



The function t(s) could be obtained as 



t(s) = -2(1 + d + e)s + (1 + 2d), d= y/e 2 - /3 2 , 



(60) 



if 7r(s) = — (d + e)s + d is chosen for fc_ = — (d + e) 2 . We can also find the eigenvalues from 
Eqs. (6) and (10) as 



A = - (d + e) (d + e + 1), A n = 2n (d + e + 1) + n(n - 1). 



Thus, solving Eq. (61) for the energy eigenvalues, we obtain 



(61) 



E nq (V ,ia) 



H 2 



2ma? 



n + l 



7 



(n + l) 



ma 2 V 

, 7 = — § — ? < n < oo. 
h q 



(62) 



On the other hand, the radial wavefunctions in the present case become 



ip ng (s) = N nq s d (l - S yPi 2d ^(l - 2s), (63) 
with s(r) = (1 + qe -Hr-R°)/a}-i and ^ 

is a new normalization constant determine by 
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= 2 (n + 2e)\T(n + 2d+lf f " (-l)y-* 

' "" l " T(n + 2rf + 2e + 1) p!(ra - p)!(2e + p)\T(n + 2d - p + 1) 



x 



rl(n-r)ir(2d + r + l) J ^ r) ' (64) 



l 

where the integral I nq (p,r) = J s n + 2d + r -p(\ _ qs) p+2e ds is given by 

o 



4 9 (P,7") =2 Fi(n + 2c + r-p+l,-p-2e : n + 2c + r-p + 2 : l)B(n + 2c + r -p + 1, 1), (65) 

Figures 1 and 2 show the variation of the ground-state (i.e., n = 0) as a function of the 
coupling constant Vo for different positive and negative q, and a = A c . Obviously, in Figure 
1, the non-Hermitian PT-symmetric generalized WS potential generates real and negative 
bound-states for q > 0, it generates real and positive bound-states for the same value of a 
when q < (Figure 2). Further, Figures 3 and 4 show the variation of the first three energy 
eigenstates as a function of a for (a) q = 1.0 and (b) q = —1.0 with Vo = 2.5m. Obviously, 
for the given Vo, as seen from Figures 3 and 4 all possible eigenstates have negative (positive) 
eigenenergies if the parameter q is positive (negative). It is almost notable that there are 
some crossing points of the relativistic energy eigenvalues for some Vq values. 



B. Non-Hermitian non-VT— symmetric generalized Woods-Saxon potential 

In this part, we consider two parameters {Vo, q} to be complex parameters (i.e., Vo — > iV , 
Consequently, the potential in Eq. (12) transforms to the following form 



V q (x) = Vb- p— — } - , , (66) 

1 + q 2 2 cosh (ax) — sinh(2o;a;) — 1 

which is a non-PT— symmetric but non-Hermitian. The complex energy eigenvalues of Eq. 

(66) are given by 



E, 



V 

nq ^ ± i (iV + q(n + 1)/ a) 
Iq 



\\4q 2 V 2 + {iV + q{n + l)/a) 2) ^ 
On the other hand, the corresponding radial wave functions u nq (s) for the s-wave could be 
determined as 
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u 



nq 



(s) = N nq s\l- S yPi 2b ^\l-2s), (68) 



and the lower spinor component w nq (s) is given by Eq. (42) with s(r) = (1 J riqe~( r ~ Ra ^ a )~ l . 

l 

The integral I nq (p,r) = J s n + 2b + r -p{\ — s) p+2e ds is given by 

o 



I nq (p,r) = 2 F 1 (n + 2b + r-p+l,-p-2e : n + 2b + r -p + 2;i)B(n + 2b + r -p + 1,1). (69) 

C. Pseudo-Hermiticity and VT— symmetric new trigonometric and periodic po- 
tential form 

Finally, in this part, when all the parameters {Vo,a,q} are complex parameters (i.e., 
Vq — > iV , a — > zoj, q — > ig), from Eq. (12) we obtain 

= - g2 + 2g J^ ) + 1 [? + sm(ax) + < «»(«*)] = - x). (70) 

Note that the transformed potential in Eq. (70) has a trigonometric and periodic form. 
As we remarked earlier, the above potential forms do not have any of the WS potential 
properties. Nonetheless, they have periodic, trigonometric and VT — symmetric behaviors. 
The potential in Eq. (70) is a pseudo-Hermitian potential [37,38] having a ir/a phase 
difference with respect to the potential (I), it is also a VT — symmetric, rj = P-pseudo- 
Hermitian (i.e., PTV^x^PT)- 1 = V q (x), with P = rj : x -> £ - x and T : % -> -i) but 
non-Hermitian having real spectrum given by 



E nq = -^ T (V + q(n + l)/a) 



m 2 



\2 ' 



(71) 



25 V v 2 -(\/ + g(n+l)/a) ; 

if the same restrictions after Eq. (55) are achieved. 

On the other hand, the corresponding radial wave functions u nq (s) for the s-wave could 
be determined as 

lM«) = N nq s c (l - sfPf^)(l - 2a), (72) 

and the lower spinor component w nq (s) is given by Eq. (48) with s(r) = (1 +iqe~ l ( r ~ Ro ^ a )~ 1 . 

i 

The integral I nq (p,r) = J s n + 2c + r -p(i _ s ) p+2te ds is given by 

o 
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Inq(p,r) =2 F 1 (n + 2c + r-p + l,-p-2ie : n + 2c + r-p + 2;i)B(n + 2c + r-p+l,l). (73) 



IV. THE SOLUTION OF THE GENERALIZED WS POTENTIAL FOR q = 

We have obtained the bound state solutions of the generalized WS potential with q ^ 
and the explicit form of the eigenvalues and the spinor wavef unctions. In addition, we start 
finding solutions for q = case with the definition of new variable s = e~ ax . Hence, Eq. 
(18) is reduced to the generalized equation of hypergeometric type: 

d 2 u n (s) 1 du n (s) 1 



+ + 



ds 2 s ds s 2 
Further, we define the following dimensionless expressions: 



e 2 + (3s + 7s 2 u n (s) = 0. (74) 



e = -(E -m), fi = i— + — — , l=—o, 75 



f(s) = 1, cr(s) = s, or(s) = -e 2 + /9s + 7s 2 , (76) 

with real e 2 > (-E 2 < m 2 ) for bound states. The substitution of the above expressions into 
Eq. (9), together with a'(s) = 1, gives 



7r(s) = ±i\J^s 2 — (k — f3)s — e 2 . 
Substituting the value for each into the last equation, we obtain 



tt(s) = ± < 



ids + e; for k = f3 + 2?5e, 
i5s — e; for k = /3 — 2i8e, 

where 5 = yfy = ^. Therefore, the most suitable form is established by 



(77) 



(78) 



7r(s) = — ids + e, k = (3 — 2i5e. 



(79) 



The following track in this selection is to achieve the condition r(s) < in Eq. (5), which 
can be obtained as 



r(s) = -2i5s + 1 + 2e, r(s) = -2i<J, \ = @ — 2i5e - iS. 
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(80) 



Let us calculate the wavefunctions. Inserting, n(s) and a(s) in Eq. (4) and consequently 
solving the resulting first-order differential equation, we find 



(s) = s e e- iVos/a . 



(81) 



In addition, to find the function, y nq (s), which is the polynomial solution of hypergeometric- 
type equation, we multiply Eq. (3) by p{s) so that it can be written in self-adjoint form 
[28] 



sy'^ + [1 + 2e - us] y' n 



y n = 0, v = 2iS. 



(82) 



The solution can be written in terms of confluent hypergeometric function as follows: 



n ( 1 -P<* , « 2iV " 

y n (s) =i F 1 I e + - + z— , 1 + 2e, 



(83) 



and consequently the upper spinor 



w(s) = A lFl (e + I + igL, 1 + 2c, *«e-«W°. 



(84) 



Finally, the lower spinor is found by 



mw(s) 



i— + E + V s 

ax 



2V 



a 



(85) 



V. RESULTS AND CONCLUSIONS 



In this work, we have seen that the s-wave Dirac equation with the generalized WS 
potential can be solved exactly for its bound states using the hypergeometric method. 
The relativistic bound states energy spectrum and the corresponding wave functions for 
the generalized WS potential have been obtained by the hypergeometric method. Some 
interesting results including the VT— symmetric, non-PT— symmetric non-Hermitian, and 
non-PT— symmetric P-pseudo-Hermitian versions of the generalized WS potential have also 
been discussed for bound states. In addition, we have discussed the relation between the 
non-relativistic and relativistic solutions and the possibility of existence of bound states for 
complex parameters. 
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FIG. 1: The ground-state (n = 0) energy, in a non-Hermitian PT-symmetric potential given by 
Eq. (54), as a function of the coupling constant Vq for three positive shape parameters q with 
a = A c . 

FIG. 2: The ground-state (n = 0) energy, in a non-Hermitian PT-symmetric potential given by 
Eq. (54), as a function of the coupling constant Vq for three negative shape parameters q with 
a = A c . 

FIG. 4: The first three energy eigenstates, in a non-Hermitian PT-symmetric potential given by 
Eq. (54), as a function of the range parameter a for a positive shape parameter (q = 1.0) with a 
coupling constant Vq = 2.5m. 



FIG. 3: The first three energy eigenstates, in a non-Hermitian PT-symmetric potential given by 
Eq. (54), as a function of the range parameter a for a negative shape parameter (q = —1.0) with 
a coupling constant Vq = 2.5m. 
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